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1—5 ; 

l^«^ ' We establish weak convergence of the empirical process on the 

*vj ■ spherical harmonics of a Gaussian random field in the presence of an 

' unknown angular power spectrum. This result suggests various Gaus- 

1 1 1 sianity tests with an asymptotic justification. The issue of testing for 
P^ ■ Gaussianity on isotropic spherical random fields has recently received 
^0 ' strong empirical attention in the cosmological literature, in connec- 

(-H , tion with the statistical analysis of cosmic microwave background 

■ 4-^ ' radiation. 

C^ ■ 

g: 

i__i, 1. Introduction. In recent years an enormous amount of attention has 

been devoted to testing for Gaussianity for splierical random fields, espe- 

^ ■ cially in the astrophysical and cosmological literature. The empirical mo- 

\Q ', tivation for these studies can be explained as follows. The ongoing NASA 

^N ■ satellite mission MAP and the forthcoming ESA mission Planck will probe 

vQ . cosmic microwave background radiation (CMB) to an unprecedented accu- 

C^ ', racy. CMB can be viewed as a signature of the distribution of matter and 

^ ■ radiation in the very early universe, and as such it is expected to yield very 

^^ . tight constraints on physical models for the Big Bang. For the density fluc- 

1 -^ I tuations of this field, the highly popular inflationary scenario [see Peebles 

C^ ' (1993) or Peacock (1999)] predicts a Gaussian distribution, whereas alter- 

f^ . native cosmological theories, such as topological defects or nonstandard in- 

^ I flationary models, predict other types of behavior. The density distributions 

of fluctuations are also instrumental for drawing correct inferences on the 

physical constants which can be estimated from CMB radiation; indeed, 

C^ I point and interval estimation procedures for cosmological parameters have 

been based almost exclusively upon Gaussian assumptions, which of course 
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2 D. MARINUCCI AND M. PICCIONI 

need to be validated before reliable statistical inference can take place. For 
these reasons, many different Gaussianity tests were considered in the recent 
cosmological literature, some of them based upon the topological properties 
of Gaussian fields [Novikov, Schmalzing and Mukhanov (2000), Phillips and 
Kogut (2001) and Dore, Colombi and Bouchet (2003)], others on higher- 
order cumulant spectra [Winitzki and Wu (2000) and Komatsu and Spergel 
(2001)]. 

Let T{6,ip) be a random field indexed by the unit sphere 5^; that is, for 
each azimuth < ^ < vr and elongation < (p < 27r, T{9, if) is a random vari- 
able defined on some probability space. Throughout this paper we assume 
that T{0,ip) has zero mean, finite variance, is mean square continuous, and 
is isotropic, that is, its covariance is invariant with respect to the group of 
rotations. Furthermore, we introduce the spherical harmonics, defined by 



21 + 1 (/ — 77i)' 

YUO,^) = \i 4vr ^i^^y PU^o.e)eM^m^l for m > 0, 

{-irYi%,XO,v), form<0, 

where the asterisk denotes complex conjugation and Pira{cos9) denotes the 
associated Legendre polynomial of degree /, m, that is, 

^Ka^) = ^^(^'-1)'> m = 0,1, 2,...,/,/ = 1,2,3,.... 

A detailed discussion of the properties of the spherical harmonics can be 
found, for instance, in Liboff (1998), Chapter 9, and in Varshalovich, Moskalev 
and Khersonskii (1988), Chapter 5. The following spectral representation 
holds in the mean square sense [see, e.g., Hannan (1970), Wong (1971) and 
Leonenko (1999)]: 

oo I 

(1) Tie,ip) = J2 E aimYimiO,V>), 

1=1 m=~l 

where the triangular array {aim} represents a set of random coefficients 
which can be obtained from T[0, ip) through the inversion formula 

(2) ai^= r rT{e,ip)Yi*r^{e, if) sine dOd^, 

Jo JO 

m = 0,±l,...,±/,/ = l,2,.... 

These coefficients are zero-mean and uncorrelated [see Wong (1971), pages 
253 and 254]; hence, if T{9, ip) is Gaussian, as we shall always assume in this 
paper, they have a complex Gaussian distribution, and they are independent 
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over I and ttt, > (although ai^^m = (~l)™0;*m)' with variance E\aim\'^ = Ci, 
m = 0, ±1, . . . , ±1. The sequence {Ci} denotes the angular power spectrum 
of the random field; we shall always assume that Ci is strictly positive for 
all values of /. 

Our purpose in this paper is to construct Gaussianity tests based on the 
empirical distribution function for the triangular array of random coeffi- 
cients {aim}, I = I, ■ ■ ■ ,L, m = 0, ±1, . . . , ±1, as L —> oo. In principle, expres- 
sion (2) can be used to recover the coefficients aim for any value of /. In 
practical applications, however, L is finite and determined by the resolution 
of the experiment. For instance, the coefficients aim niay be contaminated 
by noise, with a decreasing signal to noise ratio as I grows. Throughout this 
paper we adopt the simplifying assumption that the noise is negligible only 
for / = 1, . . . , L, hence the higher-order coefficients are discarded. The cur- 
rent status of technology in satellite missions suggests that our assumption 
may provide a reasonable approximation for L as large as several hundreds 
for the ongoing NASA experiment MAP, a value which is likely to increase 
to a few thousands for the forthcoming ESA experiment Planck. 

Testing Gaussianity in a spherical field poses considerable extra difficulty 
vs. testing in Euclidean spaces, since (fixed radius) spherical fields are never 
ergodic. Our proposal is to consider the asymptotic behavior of the empiri- 
cal process on the triangular array {a/m}) in the presence of infinitely many 
unknown parameters {Ci}\ in fact, the angular power spectrum is typically 
unknown in practice, and needs to be estimated nonparametrically from 
the data [Miller, Nichol, Genovese and Wasserman (2002) and Wasserman, 
Miller, Nichol, Genovese, Jang, Connolly, Moore, Schneider and the PICA 
Group (2001)]. The presence of a growing sequence of estimated parameters 
makes it hard to exploit the modern theory of empirical processes, as pre- 
sented for instance by van der Vaart and Wellner (1996) or Dudley (1999); 
we hence resort to a more traditional approach, based upon standard weak 
convergence theorems in the Skorohod space Z)[0,1]^, as presented for in- 
stance by Bickel and Wichura (1971) or Shorack and Wellner (1986). We 
refer to these works for the definition and topological properties of such a 
space. For the sequel it is enough to recall that, in order to prove the weak 
convergence of a sequence {Ki{r, a)} of L'[0, l]^-valued processes to the field 
K{r,a), we need to prove both convergence of all finite-dimensional distri- 
butions and tightness. For the latter we will use a sufficient criterion due to 
Bickel and Wichura (1971). 

First, for a generic "block" B = (ai,a2] x {fi,r2\ C [0, 1]^, define the in- 
crements of the field Ki, 

Kl{B) = KL{a2,r2) - i^L ("2 , y'l ) - KL{ai,r2) +KL{ai,ri). 

We define two types of adjacent blocks: 
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Type I blocks 

Bi = (ai,a] X (ri,r2], 

B2 = (a, 02] X {ri,r2], < ai < a < 02 < 1,0 < ri < r2 < 1, 

and Type II blocks 

Bi = (01,02] X (ri,r], 

B2 = (01,02] X (r,r2], 0< oi <02 < l,0<ri <r <r2 < 1. 

Bickel and Wichura (1971, Theorem 3) show that tightness of the se- 
quence {Kl} is satisfied if there exist /3 > 1, 7 > such that for all blocks 
Bi and B2, 

E{[mm{\KL{Bi)\,\KL{B2)\}r) < C{f,{B^ U ^2))^ 
which is implied by the stronger condition 

(3) E{\KL{Bi)\'''\KL(B2W')<C{fi{BiUB2)f, 71 + 72 = 7, 

where fi is some finite measure on [0, 1]^ with continuous marginals. They 
also show that for some particular class of processes (partial-sum processes) 
we can restrict to blocks having corners in [0,1] x {0, j^, . . . , 1}, so that 
(r2 — r),{ri — r)> -^ always [see Bickel and Wichura (1971), page 1665]. 

The plan of this paper is as follows. In Section 2 we start from the case 
where the angular spectrum is known and then provide a formal definition of 
the empirical process for spherical harmonics we are interested in, for which 
we state the main weak convergence result. Section 3 presents the main steps 
of the proof; Section 4 draws some conclusions, which are also illustrated by 
a small Monte Carlo experiment, and points out some directions for further 
research. Many technical lemmas are collected separately in the Appendix. 
In the sequel, we use C to denote a generic, positive and finite constant, 
whose value may vary from line to line. Also, throughout the paper we 
define — log x = 00 and (by continuity) x log x = x log x = for x = 0. 

2. The empirical process with unknown angular power spectrum. We 

start by assuming that the sequence of coefficients {C/}/=i^2,... in the angular 
power spectrum of T{9,ip) is known. Now recall that, under Gaussianity, 
|a;oP/C/ and {210/^1^}/^ = {'^\0'i,~m.\'^}/Ci are mutually independent chi- 
square variables, with one and two degrees of freedom, respectively. The 
special distributional properties of the single term |a/oP greatly complicates 
notation, whereas it can be shown that this term has no effect on asymptotic 
distributions; hence in the sequel we shall simply drop it and focus on |a/mP 
for m= 1,2,...,/. First introduce the Smirnov transformation 

/ Ig |2\ 

(4) ui„i = l-exp( ?^)' m=l,2,... ,1,1 = 1,2,. ..,L, 
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to convert the random variables aim. to a triangular array {uim} of i.i.d. 
random variables with a uniform distribution in [0,1]. We can hence define 
their empirical distribution function over the ^th row, 



1 ' 



Fi{a) = -Y,Huim<a), 0<a<l 



m=l 



KL{r,a) 



[!(•) denoting the indicator function], and the empirical process 

Gi{a) = Vl{Fi{a) - «} = ij I] l(n/™ <a)-a\, < a < 1. 

In order to detect departures from the Gaussianity assumptions over some 
region of the angular decomposition, we shall consider the doubly indexed, 
integrated empirical process 

[Lr] 

Gi{a) is zero-mean with covariance function (ai A a2)(l — (ai V a2))i as is 
well known. The integrated process -fCi(r, a) has independent increments 
in r, from which it is easily obtained that its limiting covariance function 
should be 

lim EKL{ri,ai)KL{r2,a2) = (ri Ar2)(ai Aa2)(l - (ai Va2))- 

L— >oo 

The Gaussian zero-mean process K with this covariance function is called 
the Kiefer-Miiller process on [0, 1]^. In fact it is a standard result that Ki, 
whose sample paths clearly belong to -D[0, 1]^, converges weakly to this field 
as L ^ oo. We will be concerned, however, with the much more interesting 
(and difficult) case when the power spectrum is unknown and we have to 
estimate it nonparametrically from the data. A natural candidate to replace 
Ci is the maximum likelihood estimate 



^ 1 ' 



ra=l 

The uim are then replaced by the random variables 

/ la P\ 

y«m = 1 - expf ^^ ) ^ ^ ~ exp(-/^;m), 

V Ci / 

\n |2 
t _ \"'lm\ ™ _ 1 9 / 

Qm — ; I 12 ) '"' — J^) -^1 • • • ) '■; 

where {^n, . . . ,^ii) no longer has independent components, but has a Dirich- 
let distribution with parameters (1, . . . , 1), that is, it is uniformly distributed 
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on the unit simplex of R . For the sequel we recall that the ^u, i = 1,. . . ,1, 
are exchangeable with marginal distributions given by 

P{^ii<a} = l-{l-ay-\ 



(5) 
and 



ae[0,l], 



(6) Ptel<«l,62<«2} 


= i_(i_«^y-i_(i 


-a2)'- 


' + 


^-ai-a2)'-^ 






Qi,a2G[0,l]. 


We define now 






(7) 

1 


> ; {Hvim <a)-a} 

.m=l 
_m=l ^ ^ 


-a)\ 


-.}] 


5 



a representation which will be useful for the arguments to follow. Notice 
that, as usual, G;(0) = Gi{l) = 0. The process (7) is known in the literature 
as the normalized uniform spacings process; see, for instance, Shorack and 
Wellner (1986), pages 731-733, where its limiting behavior as / ^ oo is de- 
rived. More precisely, it is shown that Gi converges weakly in D[0, 1] to Goo, 
where Goo is a mean-zero Gaussian process with covariance function 



(8) 



£'Goo(ai)Goo(a2) = ("i A a2){l - ("i V 02)} 



- (1 -ai)(l -a2)log(l -ai)log(l-a2). 
We shall focus instead on the partial sum of empirical processes 

^ [Lr] 

KL{a,r) = —=y^Gi{a), 0<a<l,0<r<l. 



As we shall see, the asymptotic behavior of K^ as L ^ 00 does not follow 
trivially from the asymptotics for Gf, intuitively, this is due to the effect of 
higher-order terms, which have to be controlled as they are summed over /. In 
particular we shall now show that the process K^ has a nonnull asymptotic 
bias. Define 

log(l - a) 



so that 



bi{a) = lEUUn< 
EGi{a) = 



Vl 



I 



bi{a). 



a 
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Lemma 2.1. As I ^oo, 
(9) lim bi{a) = (1 - a) log(l - a) + i(l - a) log^(l - a) = b{a) 



and also, as L ^ oo, 



lim EKL{a,r) = 2^/^b{a). 

L—*oo 



Proof. From (5) and for I suitably large, 

Therefore 

, . X ,f / log(l-a)\'-i 
bi{a) = lh-(l + ^^^ '-] -a 

Using Lemma A. 2, we have immediately (9). Also 

1 ^"^''' 1 h( ) ^"^''' 1 



Now 



lin^'H^l 






l^oo i/^ f^iVl Jo \pU- 

Also, because lim^^oo 6;(q;) = 6(a), for any (^ > there exists /q such that, 
for all I > Iq, \bi{a) — 6(a) | < S; hence 



-Lf2^{biia)-b{a)}<\ sup \b,ia)\ + \b{a)\]^+2SV^<e, 



OHaj| + |o^aj|f> 

for any e > 0, for L large enough, because b is arbitrary. Thus the proof is 
completed. D 



We shall hereafter write for brevity 

log(l — a) , . , 

ai = ^^ for a G [0,1]. 

The limiting behavior of the covariances was anticipated in (8); indeed, we 
have the following. 
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Lemma 2.2. For all < ai, 02 < 1, 
lim CoY{Gi{ai),Gi{a2)} 

(10) =(ai Aa2){l-(ai Va2)} 

- (1 - ai)(l - 02) log(l - ai) log(l - 02), 
and as L ^ co, /or a// < ai, 02 < 1, < ri, r2 < 1, 
lim Cov{KL(ai,ri),KL(a2,r2)} 

L— >oo 

(11) = (ri Ar2)[(ai Aa2){l - (ai Va2)} 

- (l-ai)(l-a2)log(l-ai)log(l -02)]. 

Proof. The first limiting result is standard in the theory of spacings 
[see Shorack and Wellner (1986) and Shorack (1972)]; indeed, from (5), (6) 
and Lemma A.l the reader can check that 

lim Cov{l{Cii <aii), 1(^/1 <a2i)} 

t— >oo 

= (ai A a2){l - (ai V 02)}, 
lim /Cov{l(6i <aii), l{^i2 <a2i)} 

t— >oo 

= -log(l -ai)log(l -a2)(l -ai)(l-a2), 

whence (10) follows easily. On the other hand, (11) is an immediate conse- 
quence of the Kronecker lemma and (10). □ 

Now write 

Kl{a,r) = KL{a,r)-2^bia), 

and define K*(a,r) as the zero- mean Gaussian process on [0, 1] x [0, 1] with 
covariance 

EK*{ai,ri)K*{a2,r2) 

= (n A r2)[(ai A a2){l - (ai V 02)} 

- (1 -ai)(l -a2)log(l -ai)log(l -02)]. 

To the best of our knowledge, the field K*{a,r) has not appeared in the 
literature so far, and we label it a modified Kiefer-Miiller process. Our main 
result will be the following. 

Theorem 2.1. As L^ 00, weakly in D[0,1]'^, 

=^ denoting weak convergence in the Skorohod space D[0,1]^. 
The proof of this result will be given in the next section. 
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3. The weak convergence proof. 

Proposition 3.1. AsL^co, 

Proof. Consider < ai < • • • < a^ < 1, < ri < • • • < r^ < 1, and let 

(12) Gt{a,/3]=Gi{a,l3]-EGi{a,p], 

for < a < /3 < 1. Because of Lemma 2.1 the convergence of the s x t vector 
with components K^{ai,rj) with i = 1, 2, . . . , s, j = 1, 2, . . . , i, to the same 
components of K*{a, r) is equivalent to the joint convergence of the centered 
increments 

KL{iai,ai+i] X {rj,rj+i]) - EKL{{ai,ai+i] x {rj,rj+i]) 

-, [Ln+i] 
= 7f E [Gi{ai+i)-Gi{ai)-E{Gi{a,+i)-Gi{a.i)}] 

^^ l=lLn]+i 

for i = 0, . . . , s and j = 0, . . . , t, where we have set ao = tq = 0, a^+i = Tt+i = 
1, to the increments 

K*((ai,ai+i] X (rj,rj+i]) 

= /i:*(ai+i,rj+i) + /s:*(ai,rj) - K* {ai,r ^j^^ - K*{ai+i,rj). 

Because of the independence over / we can restrict ourselves to a fixed in- 
terval for r, by simplicity (0,1], say. It is now clear that we have to use 
the multidimensional central limit theorem (CLT) with independent but not 
identically distributed summands. A suitable control on the fourth moments 
clearly implies the Lindeberg condition, namely 

1 ^ 

for any i,k = 0, . . . ,t, possibly equal. Now in Proposition 3.2 it will be shown 
that EG^ {ai, ai+i]'^ G^ {ak, ak+i]"^ is uniformly bounded by a constant, which 
will allow us to complete the proof. D 

Proposition 3.2. For L = 1,2,... , the sequence of fields Kl{a,r) is 
tight in Z) [0,1] 2. 
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Proof. We write 

Kl{a,r) = KL{a,r) - EKL{a,r) + KL{a,r) 

-KL{a,r) + EKL{a,r) -2^/¥b{a), 
where 

(14) " ' ~ 



KL{a,r):=^=J2Giiti{a)), 



, / N i a, for a<l-l ^/^, 
*'^") = li; fora>l-r3/2: 

Notice that the modification Gi{ti{a)) differs from the integrated empirical 
process Gi{a) because it is tied down to zero at a = 1 — 1~^'^, I = 1,2, ... ,L. 
The result will be clearly established if we can prove that the sequence 
KL{a,r) — EKL{a,r) is tight, and as L — > oo, 

(15) s\xg\KL{a,r) - KL{a,r)\ =Op(l), 

(16) SMp\EKL{a,r)-2y/^b{a)\=o{l). 

The proofs of (15) and (16) are given in Lemmas A. 4 and A. 5, respectively. 
We shall hence focus on tightness; let us write 

K'L{a,r) =KL{a,r) - EKL{a,r). 

It is sufficient to establish that, for some probability measure /i(-) with 
continuous marginals, 

E{\K'^{B,)\'\K'^{B,)\') < C{^^{B, U B,)f, 

where Bi,B2 are either Type I or Type II blocks. Let us consider Type II 
blocks first, for which 

E{KLi{ai,a2] x (ri,r])^/s:[((ai,a2] x {r,r2]f} 

(17) _ _ 

= E{K'L{{aua2] x (ri,r])'}i?{Ki((ai,a2] x {r,r2]f}. 

Now recalling (12), 



E{K'L{{ai,a2]x{ri,r]f} = j ^ EGf{{ai,a2]n{0,l - r'/']). 
For 0<a </?<!, 



^ l=[Lri]+l 



1 ' 






m=l 
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where Z/„,(a,/3] = ZimiP) - Zim{a) and 



/ J V I 

Thus 

0, for ai > 1 - /-^/^ 



.BGf ((ai,a2A(l-/-3/2)])^ otherwise. 
For < a < 1 write Ti{a) = a A (1 - T^/^), so that 
EGf{ai,Ti{a2)] 

= T 51 X! -E^/mi(ai,Ti(a2)]^/m2(aii'^/(a2)] 
mi=l 7712=1 

= SZfi(ai,ri(a2)] + {I - '^)EZii{ai,Ti{a2)]Zi2{ai,Ti{a2)]. 
Now note that 

EZf^{a,Ti{a2)] 

= Var(]l(ai < y^^ < Ti{a2))) < El{ai < yim < Ti{a2)) 
= Pi{ai,Ti{a2)) <e^|r;(a2) -ai| <C\a2 -ai|, 
the last step fohowing from Lemma A. 3. Also, from Lemma A. 6, 

{I- l)EZii{ai,Ti{a2)]Zi2{ai,Ti{a2)] <Cq{ai,Ti{a2)) <Cq{ai,a2) 
with 



/•02 

q{ai,a2)= (1 + | log(l - y)|) dy; 



indeed, 

q^{a,b)<Cq{a,b) for all < a < 6 < 1, 
where C does not depend on a or b. Hence 

E{^i((ai,a2]x(ri,r]f} 

<— ^ {|a2 -ail +g(ai,a2)} 

i=[Lri]+l 

= C {\a2 -ail +q{ai,a2)]. 
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Note that, since we can restrict to (r2 — ri) > -j^, 

\Lr] - \Lri] 

This bound obviously holds for B2 as well, hence 

i?(K£((ai,a2] X {n,r\fK'L{{ai,a2\ x (r,r2])') 
< C(r2 -ri){\a2 -ai| + g(ai,a2)}^- 
Let us now consider Type I blocks: we have 

E{K'L{{ai,a] x (ri,r2])'i?i((a,a2] x [n^r^]?] 

-I [Lr2\ 

= 72 E i^Gr2((ai,a]n(0,l-r3/2]) 

i=[Lri]+l 



xGf2((a,a2]n(0,l-r3/2]) 

+ 7 E i5;Gf((ai,a]n(o,i-r3/2]) 

I «=[Lri]+l J 

X 7 E i?Gr'((«,«2]n(o,i-r3/2]) 

I ^ i=[Lri]+l J 

+ - Y. i^G'r((«i,a]n(o,i-r3/2]) 

I «=[Lri]+l 



xGr((a,a2]n(0,l-r3/2])l . 



The middle term is the same as for Type II blocks, while the last term is 
bounded by the former because of the Cauchy-Schwarz inequality. As far as 
the first term is concerned, notice that 

EGf{{aua\ n (0, 1 - r^i'^])Gf{{a, 02] n (0, 1 - r^'"^]) = 0, 

if a > 1 — l^^'"^ . Otherwise, the right-hand side is equal to 

EGf{ai,a]Gf{a,Ti{a2)] 

(18) =-^EZl{ai,a]Zl{a,Ti{(^2)] 

(19) +2——EZ'i^{ai,Q\Zii{a,Ti{a2)]Zi2{a,Ti{a2)] 

(20) +2——EZii{ai,Q\Zf^{a,Ti{a2)\Zi2{ai,Q\ 
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(21) + ^{l -2)EZf,{ai,a]Zi2{a,Tiia2)]Zi3{a,Ti{a2)] 

(22) + ^{l -2)EZf^{a,Tiia2)]Zi2{aua]Zi3{ai,a] 

(23) +4:——EZii{ai,a]Zii{a,Ti{a2)]Zi2{ai,a]Zi2{a,Ti{a2)] 

(24) +Lzli?z2(ai,a]z2(a,rKa2)] 



(25) 



(26) 



+ 4^(^-2) 

X EZii{ai,a\Zii{a,Ti{a2)]Zi2iai,a\Zi2{a,Ti{a2)] 
+ ^(/-2)(/-3) 



X EZii{ai,a]Zi2{ai,a]Zi3{a,Ti{a2)]Zi4{a,Ti{a2)]. 
Now for (18) note that 

l(ai < yii < a)t{a < yn < 02) = 0, 
(27) 

pF(a,/3)<pK«,/3)<l, 
and then 

EZf^ {ai,a\Zf^ (a, r/(a2)] 

= E{l{ai <yii <a) -pi{ai,a)}'^ 

X {l{a<yi2 <Ti{a2)) - Pi{a,Ti{a2))}^ 

< Cpi{ai,a)pi{a,Ti{a2)) < C\a - ai\ \a2 - a\. 

For (19), in view of (27), simple manipulations and Lemma A. 6, 

E{{t{ai < yn < a) -pi{ai,a)) Zii{a,Ti{a2)]Zi2ia,Ti{a2)]} 

= (1 - 2pi{ai,a))E{l{ai < yn < a)Zii{a,Ti{a2)]Zi2{a,Ti{a2)]} 

+ pi{ai, af E{Zn{a,Ti{a2)]Zi2ia,Ti{a2)]} 

= -(1 - 2pi{ai,a))pi{a,a2)E{l{ai <yn < a) Z 12 {a, n (02)]} 

+ pi{ai,afE{Zn{a,Ti{a2)]Zi2ia,Ti{a2)]} 

< Cpi{ai,a)q{a,Ti{a2)) <C\a — ai\q{a,a2)- 

The argument for (20), (23) and (24) is entirely analogous. For (21) we 
obtain 

EZn{ai,a\Zi2{a,Ti{a2)]Zi3{a,Ti{a2)] 
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= {l-pi{ai,a))E{t{ai < y/i < a)Zii{a,Ti{a2)]Zi2{a,Ti{a2)]} 
-pi{ai,a)E{Zii{ai,a]Zii{a,Ti{a2)]Zi2{a,Ti{a2)]} 

= (1 - 2pi{ai,a))E{Zii{ai,a]Zii{a,Ti{a2)]Zi2{a,Ti{a2)]} 
+ {1 - Pi{ai,a))pi{ai,a)E{Zii{a,Ti{a2)]Zi2{a,Ti{a2)]} 

C C 

< -rPl{oii,a)q{ai,a)q{a,a2) < y |a — ai\q{ai,a)q{a,a2), 

in view of previous results and Lemma A. 7. The argument for (22) and (25) 
is analogous. Finally, it is clear that for (26) it is sufficient to prove 

EZii{ai,a\Zi2{ai,a\Zi^{a,Ti{a2)]Zi^{a,Ti{a2)] 

C C 

< ■j^q'^{ai,a)q'^{a,Ti{a2)) < -^q{ai,a)q{a,a2), 

a result established in Lemma A. 7. Hence each of the terms in (18)-(26) is 
bounded uniformly by 

C{\a2-ai\ +g(ai,a2)}^. 

The proof can then be completed by routine manipulations. D 

Remark 3.1. It is immediately seen that the bounds of (18)-(26) given 
previously hold for arbitrary nonoverlapping blocks. As such they can be 
used to establish the Lindeberg condition (13) in the proof of Proposition 3.1. 

4. Comments and conclusions. Theorem 2.1 is immediately applicable to 
statistical inference procedures, and in particular, to tests for Gaussianity. 
For instance, a Kolmogorov-Smirnov type of test is implemented, for any 
suitably large L, if we evaluate 

(28) SL = sup\Kl{a,r)\, 

a,r 

and compare the observed value with the desired quantile of the law of 
^oo = sup^ J. \K*{a, r)\. In principle, the latter value can be derived by Monte 
Carlo simulation, since the limiting distribution does not entail any unknown 
nuisance parameter. Likewise, Cramer-Von Mises and many other types of 
goodness-of-fit statistics could be easily implemented. 

Tests for Gaussianity on spherical maps have recently been considered by 
several authors in the physics literature, as mentioned in the Introduction. 
The focus of these papers is very much on the physical discussion rather than 
on statistical methodology, so that any comparison seems inappropriate. We 
claim, however, that our method enjoys some important advantages on any 
empirical procedure in this literature. To mention a few, our procedure al- 
lows for a rigorous asymptotic theory, which is made possible by the focus on 
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harmonic coefficients; it allows for inference completely free from nuisance 
parameters, whereas in other papers test statistics are considered whose law 
depends on the values of the angular power spectrum Q. Due to our study 
of the asymptotic behavior of the whole field K^, many different testing 
procedures can be implemented; these procedures provide information not 
only on departures from Gaussianity, but also on their location in harmonic 
space; this is important, as different physical mechanisms are known to op- 
erate at the various multipoles. Moreover, the distributional properties of 
the normalized spherical harmonic coefficients may have some independent 
interest for other cosmological applications. 

It is clearly of great interest to evaluate the power of our testing proce- 
dures in non-Gaussian situations. With regard to this issue, a crucial point 
is the nature of non-Gaussianity. In the cosmological literature departures 
from Gaussianity are occasionally generated by superimposing non-Gaussian 
structures over a Gaussian map. For instance, it is possible to mimic a pop- 
ular class of topological defects models (the so-called "cosmic strings") by 
setting 

(29) T{e,ip) = T''{9,ip) + T^{e,^), 

where T [9, (p) is a Gaussian map, and T {9, (p) is a map of Poisson-distributed 
segments of randomly varying directions, length and level [for more details 
on simulations of non-Gaussian spherical fields and their physical meaning, 
see Hansen, Marinucci, Natoh and Vittorio (2002)]. T^{9,ip) and T^{9,ip) 
are taken to be zero-mean and independent; we define the percentage of 
non-Gaussianity in the map by 

Png- ^^'(^'^)' 



We fix L = 500, a realistic value for the MAP experiment, and we eval- 
uate the threshold values for sizes a = 10%, 5%, 1% by 500 Monte Carlo 
replications of (28): we obtain 0.947, 1.012, 1.160, respectively. We then gen- 
erate string maps with 100, 500 and 1000 strings, and report the rejection 
frequencies for 100 replications of model (29) (with different percentages of 
non-Gaussianity) in Table 1. 

We leave several issues for future research. To increase the power of the 
tests, it is useful to consider empirical processes constructed on several rows 
[see the applied papers by Hansen, Marinucci, Natoli and Vittorio (2002) 
and Hansen, Marinucci and Vittorio (2003)]. As these rows are independent 
under the null of Gaussianity, the extension is conceptually straightforward, 
although computationally burdensome. For cosmological applications, some 
additional difficulties may arise in practical implementation; in particular, 
it may be the case that T{9,ip) is subject to some measurement error, or 
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Table 1 
Power of the Kolmogorov-Smirnov test 



Size of the Test\PNG 


0.1 


0.2 


0.3 


0.4 


0.5 


100 cosmic strings 


1% 
5% 
10% 


28% 
38% 
42% 


83% 
86% 
88% 


98% 
99% 
100% 


100% 
100% 
100% 


100% 
100% 
100% 


500 cosmic strings 


1% 
5% 
10% 


6% 

14% 

20% 


41% 
49% 
56% 


82% 
88% 
90% 


92% 
97% 
99% 


100% 
100% 
100% 


1000 cosmic strings 


1% 
5% 
10% 


3% 
8% 
16% 


23% 
35% 
40% 


66% 

78% 
81% 


95% 
97% 
97% 


100% 
100% 
100% 



that it is only incompletely observed, or both. As mentioned before, with the 
current status of technology on satellite- and balloon-borne experiments, it 
is known that observational error is some order of magnitude smaller than 
signal, and thus can safely be ignored, for / as large as 1000 or more; we hence 
consider the asymptotic theory provided in this paper to be valuable for the 
practitioner, as L on the order of 10^ (with a total number of observed aim 
on the order of 10^/10^) seems sufficient for the asymptotic theory to yield 
applicable approximations. However, the approximation of these sampling 
distributions can possibly be improved by bootstrap methods, to take into 
account the nuisance parameters governing noise. The presence of gaps in 
the observed field poses more challenging questions, and will be addressed 
elsewhere. 

APPENDIX 

Lemma A.l. Let an be a sequence of real numbers such that lim„_»oo 'n'^o-n 
— c < oo. Then 

lim n((l-ha„)"-l)=c. 

n — >oo 

Proof. The proof is trivial, writing (1 + On)" = exp(?ilog(l + an)) and 
using the expansions, as r — > 

(30) log(l + r) =r- — + 0(r3), 

(31) e'' = l + r + 0{r^). 
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D 



The proof of the fohowing lemma was very much shortened thanks to 
the comments of one referee. 



Lemma A. 2. For any C >0, as ^ oo, 

logx 



sup 
i-c<x<i 



llx- 1 



I 



X log X ■ 



X log^ X 



Proof. We have, uniformly in I ^ < x < I, 

logx\] f f n -<\i f t logx 



Hx- 1 + 



I 



dx-exp({l- l)logn + 
llx — exp f (/ — 1 



/ 



logx log X 



2/2 



+ 



log-^l 



. . logx log X ^flog I 
L<x — exp log X ; — h U 



I 



21 



P 



:/xn-exp(-i^-i^ + 0^^°^'^ 



Ixil- 1 



I 21 

logx log X 



I 



xlog^x 
: X log X H h U 



■ X log X + 



2 
X log^ X 



21 
log^/ 



+ 



/2 

log"/ 

/2 



/ 



+ o(l), 



where we used (30) and (31). □ 



Lemma A. 3. The marginal density of yim = 1 — exp{—lS,im.) is hounded 
uniformly by e^ for l>2. 

Proof. The marginal density of S,im is 

f^,jt) = {i-i)ii-ty~H[o,i]{t), 

whence by the change of variable formula, 



Jyim \y) 



I -I ( log(l-j,) y-2 



/(1-y) 



I 



^ I (, ^ \og{l-y) \'~\ ,. 

<- -11 + : I l[0,l-cxp{-0)(2/)- 



(1-y) 



I 
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Now let X = 1 — y and differentiate with respect to x; we obtain 
(1 - 2//)(l + r^ loga;)'-3 _ (1 + rMogx)'~2 

{1 + l^^logxy-^l-^{2 + logx) 



x^ 



which is equal to zero at x = e ^; the latter is easily seen to be a unique 
maximum. Hence 






2\ «-2 

D 



Lemma A. 4. AsL^oo, 

sup\KLia,r) - KLia,r)\ =Op{l). 

a,r 

Proof. By the definition of KL{a,r) (14), 

_ ~ 1 ^ 1 ' 

sup|iCL(a,r)-Ki(a,r)|<^=^^^l(l-r3/2<y,„<l); 

hence, using Lemma A. 3 and Markov's inequality, 
sup|KL(a,r)-KL(a,r)|=Op( VyiS{l(0<e^'^'- < T^/'^)}] 



D 



In the sequel, recall that 

Zim{{a,/3]) = l{ai < yim < /3) -pi{a,/3), 
pi{a, P) = El{ai <yim<(3)- 

Lemma A. 5. AsL^oo, 

sup\EKL{a,r) - 2y/?b{a)\ =o(l). 
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Proof. By the definition of KL{a,r), 
EKL{a,r)-2^/^h{Q.) 

[Lr 



1 



^/^S 



-rMtiia)) 



2^/^b{a) 



1 K 1 



wliose absolute value is bounded by 

(32) _L^i=sup|6KtK«))-&( 

(33) + sup \b{a)\ sup 



a] 



a) 



1 K 1 



vig;?!-'^. 



1 K 1 



Now for (33) we have that sup„ I^(q^)I ^ C, whereas approximating the sum 
with the integral, it is easy to see that 



^ [Lr] ^ ^ [Lr] ^ 



VLfr[Vl Lf^^ ^/TfL ~ Jo Vx 



r I 

< I —=dx = 2^/^, 



and 



Jo 



1 



i/L I r ^ 

—= dx + —= dx 

\/x Jl/L \Jx 



2 l!^ 



1 



L L 



1=2 



Hence 



(33)<C-^ = o(1) 



V(^-i)/i 



as L — > 00. 



On the other hand, for (32) it is enough to prove that 
sup \bi{ti{a)) -b{a)\ 

0<a<l 

< sup \bi(ti{a)) — b{a)\ + sup \b{a)\ 

0<a<l-«-3/2 i_;-3/2<Q,<]^ 

is 0(1) as I -^ cxo. For the second term on the right-hand side, continuity of 
b{a) at a = 1 is enough. For the first, we take x = 1 — a, so that we need to 
establish that 



lim sup 



Ux 



1 + 



logx 



l~l 



1 



x log X — X log X 



0, 
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which fohows from Lemma A. 2. Thus the proof is completed. D 

Lemma A. 6. For < ai < 02 < I - /"^/^ < /?i < /32 < 1 - l~^^^, we 
have 

C 

I Cov{l(ai < yii < 02), l(/3i < Vn < /52)}| < jg(ai, a2)g(/?i,/32), 



where 



q{a,b)= f {l + |log(l-y)|}dy 

Ja 
rl—a 

= {1 + |logx|}(ix 

Jl~b 



n~b 
is a finite measure and C does not depend on I, a, (3. 

Proof. After a change of variable we have 

rl—ai f-l~l3i 

Cov{l{ai<yii<a2),l{(3i<yi2<P2)}= / hi{xi,X2)dxidx2, 

Jl~a2 Jl-132 

where the covariance density /i;(xi,X2) is given by 

/-3 



1 \ {l-l){l-2) ( logxiX2 Y 
1^2 1 l^ V I J 



for 



Xl3;2 

_ (/ - 1)' A ^ iogx^y-2 / ^ log 3:2 y-^ 

= Ai{xi,X2) X Bi{xi,X2), 



A I \ U-\ ^ l0gXlX2\V logXi\2/ logX2\^\ ^ 



1 /'^ ^ logXlX2\ 
X1X2 



X— |1 + ^^J 



2/ i„„^_\2 



(; - 1)^ (1 + ri iogxi)2(i + /~i iogx2) 

^2 (l + /-llogXlX2)' 

Now for Ai{xi,X2) we have 



/ logxiX2y 

(IH ^ I <2;i2;2, 
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because 








(34) 




logx flogx\ 




always. Also 








sup 


{ 


/ logXlX2\^A logxi' 

V I J \ I y 


\2/ ^logX2 


/-3/2<Xl,X2 


<1^ 




/ V I 



2~, -1 



<^ii_5MV('i_!MVr<c, 



4 / 



2 / 



uniformly for 1>1. Thus Ai (xi,X2) is bounded uniformly in x and /. Now 
call 

Cl,l — 7^5 ) C2 ; - 



P 



Thus Bi{xi,X2) becomes 



(35) 



logxi log Xi 



C2,i 1 + 2^-^ + 



1 + 2 



/ 



log3;2 , log^3;2 



/ ' P JV ' ~ I 

/ loga;ilogX2 Y 

' V l'^ + llogXiX2j 



+ 



P 



and because |c2 ; — ci /| < / ^, 



(36) 



sup 

/-3/2<a;<l 



log a; I 



<1, l^<n^ for all r^/^ < X < 1, 



/2 



^ logxilogX2 V _-^ y^ /^^\ 

/2 + /logXiX2y tiU/ 



I 

logXllogX2 \^ 
P + nogXlX2/ 



logXilogX2 Y^ //\ / logXilogX2 Y 
/ + logXiX2 ^\kJ\l'^ + llogXlX2/ 



k=2 



we obtain 



C 



1(35)1 < |c2,i -ci^i\ + — {|logxi| + |logX2| + |logxi| |logX2|} 



+ ci; 



c 



l\ f log Xi log X2 

z"^ V^y v^ + nogxiX2 



E 



< —{1 + |logXl| + |logX2| + |logXl| |logX2|} 



|logX]J|logX2| j 1 

+ ^ /-31offZ V^l 



E 

fc=2 



/\ / |logXilogX2| 
k) \P + l\ogXlX2 



k-l 
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(J 

< —{1 + |logXl| + |logX2| + |logXl| |logX2|}. 

The last step follows from 

1 ' //\ / |l0gXil0gX2| y-^ 

.^ l^^^ / I log Xi log X2 I \ ^"^ 
~i-3/4T,X2<lfr2 ^' ^2 + ^ log 2:1X2; 

_^ 1 / (9/4)log^/ \fe^^^ 
(^^k\\ Z-31ogZ ) - ' 

because log l/{l — 3 log/) is uniformly bounded in L □ 

Lemma A. 7. For all Q <ai < a2<l - r^/^, < /3i < /?2 < 1 - /"^/^ 
< 71 < 72 < 1 - /"3/2^ < 5i < ^2 < 1 - /~^/^ we /iat;e 

I EZii {ai,a2]Zi2 (/3i , /3] ^^3 (7i , 72] I 
(37) 

C 
< jg(ai, 02)9(^1, /3)g(7i, 72), 

and 

\EZii{ai,a2]Zi2{/3i,P]Zi3{-fi,-f2]Zi4{6i,62]\ 
(38) ^ 

<|2 9(ai,a2)g(/3i,/3)9(7i,72)g('5i,(52), 
where C does not depend on a, (3, 7, 5 or I. 

Proof. For the sake of brevity, we give only the proof of (38); the proof 
of (37) is similar, indeed slightly simpler, and can be found in Marinucci and 
Piccioni (2003). After a change of variable we have 

SZ/i (ai , 02] -Z';2 (/?! , ^] ^/3 (71 : 72] -Z'm (5i , (52] 

hiixi,X2,X3,X4) dxi dx2dx3dxi, 



'l—a2 Jl— /32 -'I-72 •^1-52 

where the four-term covariance density hi{xi,X2,X3,X4) is given by, for 

1-3/2 ^ xi,X2,X3,X4 < 1, 



1 


('- 


1)(/- 


m 


-3)(/- 


"'Vi 


-+■ 


logXiX2X3X4y" 
/ J 


-5 


xia;2a;3a;4 




/4 




V 




- 


1 V 


(l- 


ini- 

/4 


m- 


-3) 
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I J \ I 



'-2/ ^. . .^logXfc y-4 



x(i+i^v ri+5M>^^ 



_ 3 1 (^-1)' A + l^g^\ '" Yi + log 2^2 ^ '-' 



3:1X2X3X4 ^^ V ^ / v ^ 



«-2 / l„„^.\/-2 



x(l + i^V fl + ^^S^^ 



Define 

_ a-l)^ _ (/-lf(/-2) 

Cli ,^ 5 C21 .^ , 

(/ - If {I -2){l- 3) (Z - 1)(/ - 2){l - 3)(/ - 4) 

C3i = J4 , C4; = j^ . 

We have 

/li(xi,X2,X3,X4) 

r 4 3 4 ^ 

= 4^1 i C41C1 - Y^ C3ii]iidii + X! X! (^2i&ijidiji - "ic^iiidi \ , 

where 

tpi = Cr^ai{xiX2X3X4y~^, ai(x) = l-\ — , 

c/=(n«'(^o)(n n «f(^^^j))(n n n afix,x,xk)], 

■f]ii = n^^^^i) ) n n "F(^*^i) l n afixiXjXk) |af(a;iX2X3X4), 

%i = n '^^(^fc) n «f(^fc^i) 

Vfc^ij / \k,l^i,j J 

/ 2 3 4 \ 

^ ( II n n «f(^fca;;Xm) jaf(xiX2X3X4), 
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Cl 



and 



'34 \/234 \ 

n n "■f(^*^j))(n n n a'f{xiXjXk)]af {xix2X3X4) 

, i=ljr=i+l / \i=lj=i+lk=j+l / 



dii 
diji 



ll 



d 



P + /logXiJ;2X3X4 

{I + \ogXi){l + logXj){l + logYlk^i,j Xk) y 

/^ + ?^logj;i2;2X3X4 J 

(/ + logXi)(Z+logX2)(/ + logX3)(/ + logX4) 
l^ + P log X1X2X3X4 



Now note that 






1 + 



fc=i 

logXi{J2j^iiogXj 

/ + loga:;iX2X3X4 



logXi{J2j^i'^OgXj) 
P + / log 2:1X2X3X4 



+ eii log Xi;^ log: 



'J ' 



j¥=i 



where 



e; logx^;^logXJ■ 



E 



fc=2 



ik 



c 



l{l-l) 



' logXi{J2j^ilogXj 
P + /logXiX2X3X4_ 

logXi(J2jj^i^oSXj) Y f 27 log^ / 
l-/-i61og/ J ^^^l 4 /-61og/ 



< 



C/ 3 4 y 

72 E E |logXi||logXj| 



Also 



2/ N 1 , ^logXi log^Xi 

O; (Xi) = 1 + 2—^ + ^^ : 



logXi logXj /log^Xj + log^X 



af(x.x,) = 1 + 3^ + 3^ + 
af(x.x,x,) = l + 4i^ + 4i^ + 4i^ 

log^ Xj + log^ Xj + log^ Xfc 



^2 



+ 



/2 



,5 
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4 /v-4 - „2, 



al {xlX2X:iXi) = 1 + - 2^ log Xi + O -2 



'1=1 



where the O(-) bounds on the remainders are uniform over Xj. Counting 
terms then gives 

4 1 

= l + (2 + 3x3 + 3x4)^^|^ + p^;(xi,X2,X3,X4) 

i=l 

= 1 + 23^ 1-^+ P(l{xi,X2,X3,X4), 

C ^ 4 

|/5C/(2;i,X2,X3,X4)| < |2 XI XI |logXi| |logXj|. 



r^i; = 1 + (3 X 4 + 3 X 3 + 5) °^^'' 



+ (2 + 2x4 + 3x3 + 5) X "^^ + /0»?« (^^i; Xj ,Xk,xi) 

^OgXi ^—^ logXo 

1 + 26 — h242^ — --^ + p^i{xi;xj,Xk,xi), 



C ^ 4 

Pr,i{xi;xj,Xk,xi) < j^Y^ Y^ |logXi||log: 

i=l j=i+l 



whence 



Xr?.=4 + 98X:^ + o(^E E llogx.lllo, 



Likewise, for k,l ^ i,j, 



a 1 , na^ogxi+logxj logXfc + logx; 
6iji = 1 + 26 '- + 25 

+ <^ /2E E |logXi||logxj| , 



EE%z = 6 + 153E^ 

1=1 j=i+l i=l 

/l 3 4 



+ ^ 72E E |log^i||logXj| 
\ i=lj=i+l 
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Finally, 

6 = l + 26^i^ + ofl^ Y: llogx.lllogx.lV 

i=i '■ V* i=ij=i+i ) 

Now combining all terms, we obtain that the covariance density is bounded 
uniformly in absolute value by 

4 11 I 

I C4/ - 4C3; + 6C2i - 3cii | + | - Ic^i + 5C2i - Sci; | ^ ■ 



1=1 



4 I J 

+ |23c4/ - 98c3« + 153C2/ - 78cii| ^ ■' "^"^^ 



.=1 ^ 



+ TiE E |logXi||log2;j| . 
\' i=i i=i+i / 



Now 



C4/ - 4c3i + 6C2/ - 3ci/ 

(/ - \){l -2){l- 3){l - 4) - 4(/ - 1)2(/ - 2)(/ - 3) 



+ 



6{l-lf{l-2)-3{l-l)^ 



/4 
35/2 -4(17)/2 + 6x9-3x6 



/4 



^ + o(r3) 



Similarly, 



Iim{-2c3i + 5c2i -3ciJ 

,. -2(/ - 1)2(/ - 2)(/ - 3) + 5(/ - If {I - 2) - 3(/ - 1)4 
= hm "5 

= y + o(r2), 

lim |23c4z - 98c3i + 153c2z - 78ciz| 
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It can thus be concluded that the four-term covariance density is bounded 
uniformly by 

Cf ^ 3 4 >| 

hi{xi,X2,X3,X4)<-^ll+^\\ogXi\+Y^ Y^ \\0gXi\\l0gXj\\ 
[ i=l i=l 7=J+1 J 




1=1 \ 1=1 

and thus the proof is completed. D 
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